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Abstract.  Room  temperature  and  temperature-dependent  thermal  rate  constants  are 
calculated  for  the  state-to-state  vibrational  relaxation  of  NO^  9)  by  atomic  oxygen 
using  the  quasiclassical  trajectory  method  and  limited  ab  initio  information  on  the  two 
lowest  O+NO  potential  energy  surfaces  which  are  responsible  for  efficient 
vibrational  relaxation.  Comparisons  of  the  theoretical  results  with  the  available 
experimental  measurements  indicate  reasonable  agreement  for  the  deactivation  of 
NO(v=2,3)  at  300  K  and  NO(v=I)  at  2700  K,  although  the  calculated  relaxation  rate 
constant  for  NO(v=I)  at  300  K  is  approximately  a  factor  of  2  below  the  measured 
value.  The  state-to-state  relaxation  rate  coefficients  involve  the  formation  of  long- 
lived  collision  complexes  and  indicate  the  importance  of  multi  quantum  vibrational 
relaxation  consistent  with  statistical  behavior  in  O+NO  collisions.  The  present 
results,  combined  with  recent  measurements  of  vibrational  relaxation  for  NO(v=2,3), 
suggest  that  the  current  atmospheric  models  of  NO  cooling  rates  require  higher 
atmospheric  temperatures  and/or  an  increase  in  the  NO/O  number  densities. 
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1.  Introduction 

Fundamental  vibration-rotation  band  emission  (Av=l)  from  nitric  oxide  near 
5.3  (im  has  been  long  known  to  be  an  important  cooling  mechanism  in  the  terrestrial 
thermosphere.1-2  Vibrationally  excited  nitric  oxide  NO(v=l),  mostly  v=l,  may  be 
produced  by  the  collisions  of  ground  state  nitric  oxide  NO(v=0)  with  atomic  oxygen. 
It  may  also  be  produced,  with  vibrational  excitation  up  to  v=I8,  by  the  reaction  of 
metastable  N(2D)  atoms  with  02, 

N  (2  D)  +  02  ->  NO(v  <  18)  +  O  ( 3  P) 


which  is  exothermic  by  3.76  eV.3'5  This  reaction  may  be  the  dominant  process  for  the 
production  of  NO  in  the  natural  and  artificial  auroras.  Highly  rotationally  and 
vibrationally  excited  NO  may  also  be  produced  by  the  reaction  of  ground  state  N(4S) 
atoms  with  02, 


N  (4  S)  +  O2  NO(v  <7)  +  0  (3  P) . 


This  reaction,  although  exothermic  by  f  .38  eV,  has  an  energy  barrier  of  0.3  eV  and 
therefore  requires  translationally  hot  N(4S)  atoms  to  proceed.6  Reaction  (2)  is  an 
important  source  of  NO  in  the  quiescent  day7  and  night  time8-9  thermosphere  and 
may  also  be  important  in  natural  and  artificial  auroras.  The  fact  this  reaction 
requires  translationally  hot  atoms  to  proceed  may  also  explain  the  large  variability 
of  NO  density  with  the  geophysical  parameters.10 

The  reaction  of  N2  with  ground  state  0(3P)  atoms,  known  to  be  important  in 
combustion  where  it  produces  highly  vibrationally  and  rotationally  excited  NO,11  may 
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also  produce  NO  at  very  high  altitudes  in  the  terrestrial  thermosphere  where 
translationally  hot  O  atoms  abound,12 


N2  +  0(3  P)  -»  NO  +  N(4  S). 


The  vibration-rotation  excitation  of  NO  produced  in  the  thermosphere  by  reaction 
(3),  which  is  endothermic  by  3.25  eV,  would  depend  upon  the  excess  energy  available. 

It  is  clear  that  vibrationally  excited  NO  is  abundant  in  the  terrestrial 
thermosphere  and  the  amount  of  cooling  it  can  produce  clearly  depends  not  only 
upon  the  rate  coefficient  of  collisional  deexcitation  of  the  vibrational  level  of  the 
reactant  NO  but  also  on  the  vibrational  level  of  the  product  NO.  Vibrational 
relaxation  of  NO  due  to  collisions  with  N2  is  much  slower  than  due  to  collisions  with 
O,  therefore  only  collisions  involving  NO  and  O  have  been  considered, 

NO(v)  +  O^NO  +  Q 


where  the  vibrational  level  of  product  NO  has  not  been  identified.  The  rate 
coefficient  of  Reaction  (4)  has  been  measured  at  room  temperature  by  Fernando 
and  Smith13  and  by  Lilenfeld,14  and  for  v=2  and  3  by  Dodd  et  al.15  Surprisingly,  the 
results  of  Dodd  et  al.  for  higher  vibrational  levels  are  small  than  those  of  earlier 
workers  for  v=I  by  a  factor  of  about  2.5.  Glanzer  and  Troe,16  in  their  shock  tube 
experiment  at  an  effective  temperature  of  2700  K,  however  measure  about  the  same 
rate  coefficients  for  v= I  and  2.  A  calculation  by  Quack  and  Troe17  at  2100  K  using 
a  statistical  adiabatic  model  gives  a  lower  rate  coefficient  than  their  shock  tube 
measurement  by  a  factor  of  about  3.  However,  these  authors  assume  that  vibrational 
relaxation  occurs  only  on  the  attractive  ground  state  NQ2  potential  energy  surface.18'23 
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Recent  electronic  structure  calculations23  have  shown  that  two  surfaces  should  be 
considered  in  NO  vibrational  relaxation  by  O.  Thus,  the  previous  calculations  of 
Quack  and  Troe  should  be  considered  an  underestimate. 

The  factor  of  2.5  difference  between  the  v=f  and  v=2  rate  coefficients, 
however,  are  important  in  atmospheric  modeling  because  they  correspondingly 
impact  the  NO  and/or  O  densities.  Furthermore,  for  the  state-to-state  vibrational 
relaxation  of  NO,  there  is  only  a  single  experimental  measurement24  for  v=36vN=2 
and  one  theoretical  study  for  v=I6vN=0  and  v=26vN=0,f.17  Given  the  importance 
of  these  questions  and  inconsistencies  concerning  NO  relaxation  to  the  atmospheric 
community,  this  work  employs  a  quasiclassical  trajectory  study  to  obtain  these  rate 
coefficients  as  a  function  of  the  initial  and  final  vibrational  state  of  NO  as  a  function 
of  temperature. 

2.  Potential  Energy  Surfaces 

Although  the  spectroscopy  of  the  ground  and  excited  states  of  the  N02 
molecule  has  received  extensive  theoretical  attention  over  the  last  several  decades,18'22 
the  potential  energy  surfaces  of  0(3P)+N0(X2n)  important  for  vibrational  relaxation 
of  NO  are  less  well  characterized.  The  only  ab  initio  study  of  the  interaction  of 
0(3P)+N0(X2n)  is  the  recent  work  of  Katagiri  and  Kato,23  who  obtained  the 
asymptotic  interaction  of  the  f2  potential  energy  surfaces  (neglecting  spin-orbit 
coupling)  correlating  with  0(3P)+N0(X2n).  Including  the  effects  of  the  spin-orbit 
interaction,  18  adiabatic  doubly-degenerate  electronic  states  are  obtained.  Only  the 
two  lowest  states  are  attractive,  while  the  remaining  16  states  are  repulsive.  Since  it 
is  well  known  that  repulsive  potentials  are  characterized  by  extremely  inefficient 
vibrational  energy  transfer,  the  present  study  employs  the  lowest  two  adiabatic 
surfaces  of  I2ANand  f2AO  symmetry.  A  lack  of  any  details  concerning  the 
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interactions  among  these  surfaces  precludes  any  dynamical  investigation  of 
nonadiabatic  effects. 

In  addition  to  extensive  studies  on  the  bound  N02  molecule,  there  has  been 
considerable  work  on  the  ab  initio  2ANPES25  26  appropriate  for  N+02,  Reaction  (2), 
as  well  as  several  calculations  of  the  reactivity  using  classical  dynamics.6,25  However, 
the  ab  initio  calculations  have  focussed  on  molecular  geometries  in  the  NO  +  O 
product  channel  important  for  Reaction  (2).  The  relationship  of  the  NO+O  product 
channel  from  N+02  to  the  N02  molecule  (and  thus  O+NO  relaxation)  is  extensively 
discussed  by  Gilibert  et  al.25  where  they  have  constructed  correlation  diagrams  for 
several  molecular  geometries.  Suzzi  Valli  et  al.26  have  provided  the  most  detailed 
calculations  of  the  NO+O  product  channel,  however,  again  emphasizing  molecular 
geometries  most  appropriate  to  the  N+02  reaction.  From  these  studies,  it  appears 
that  a  barrier  of  at  least  2  eV  exists  between  the  N  +  02  reactants  and  the  N02 
molecule.  Any  evaluation  of  the  current  semiempirical  2ANPES  will  have  to  await 
more  extensive  ab  initio  calculations  at  geometries  relevant  to  Reaction  (4). 

The  analytical  form  for  the  O+NO  PES  is  obtained  by  adding  a  term  to 
describe  the  potential  wells  of  the  I2ANand  I2AO  surfaces  to  a  London-Erying- 
Polanyi-Sato  (LEPS)27  surface  (which  is  taken  to  be  the  same  for  both  electronic 
surfaces), 


V( n , r 2 > r2)  ~  V leps  ( ri , r2> r2)  +  V a ( ri , r2) 


where  rb  r2,  r3  are  the  NO,  NO,  and  02  internuclear  distances,  respectively.  The 
angular  dependence  and  long-range  behavior  of  the  potential  well  is  represented  by 
a  simple  angular-dependent  exponential  function 
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V a( ri,ri)  -  A(0 ) exp[-y{( rrr°)  +  ( ri-r)} \ 

where  r°  is  the  location  of  the  potential  well,  y  determines  the  range  of  the  potential 
well.  The  angular  dependence  is  given  by  a  3rd  order  polynomial  in  the  ONO  bond 
angle  0 


A(0)  =  ICk0k. 


The  attractive  and  repulsive  diatomic  potentials  used  in  the  LEPS  function  are 
represented  by  a  modified  Rydberg  potential28  and  the  anti-Morse  function27 


V  R(r)  =  -  De[l  +  ai(r  -  re)  +  a2(r  -  re  f  +  a3(r-re  /  J  exp{-  ai(r-  re)}  +  De 


V  AM  (r)  =  -  De  [l  +  exp{-p(r  -  re)}J  +  D, 


where  De  is  the  dissociation  energy,  re  is  the  equilibrium  internuclear  distance,  (a,,  a2, 

a3)  are  parameters  fit  to  the  RKR  potential,  and  p=a1/S/2.  In  addition  to  the 

diatomic  parameters  defined  above,  the  LEPS  surface  requires  the  specification  of 
the  adjustable  Sato  parameter,  z-„  for  each  diatomic.  Once  the  LEPS  surface  is 
determined,  the  {ck}  in  Equation  (7)  are  determined  such  that  the  angular 
dependence  of  the  potential  at  the  equilibrium  geometry  defined  by  Equations  (5)- 
(9)  reproduces  the  ab  initio  calculations  of  Hirsch  and  Buenker20  for  the  I2ANstate 
and  Gillispie  et  al.18  for  the  I2AO  state  (which  correlates  with  the  2B2  state  in  the 
strong  interaction  region23).  The  final  parameters  for  the  I2ANand  I2AOpotential 
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energy  surfaces  are  given  in  Table  1  for  the  LEPS  function  and  Table  2  for  the 
angular-dependent  exponential  function  representing  the  potential  well. 

Contour  plots  of  the  f  2ANand  f  2AOangular  dependence  of  potential  energy 
surfaces  are  shown  in  Figures  f  and  2,  respectively,  for  the  NO  bond  distance  fixed 
at  its  equilibrium  internuclear  value.  For  a  wide  range  of  bond  angles,  the  approach 
of  the  O  atom  to  the  N  end  of  the  NO  molecule  occurs  with  no  barrier.  For  the  1 2AN 
PES  the  formation  of  N02  is  favored  for  a  wide  range  of  bond  angles  (02;  85°)  for 
approach  of  the  O  atom  to  the  NO  molecule,  while  the  angular  approach  for  the  f2AO 
surface  is  more  restrictive  (f45°2;02:85o).  The  location,  depth,  and  angular 
dependence  of  the  potential  wells  for  the  bound  N02  molecule  on  both  surfaces  are 
in  qualitative  agreement  with  the  available  ab  initio  calculations.  However,  no 
attempt  was  made  to  fit  the  available  spectroscopic  data  of  the  stable  triatomic  N02 
molecule  for  either  the  ground  state  f  2ANor  excited  state  f2AO(which  correlates  with 
the  2B2  state).  The  long  range  interaction  is  in  good  agreement  with  the  ab  initio 
calculations  of  Katigiri  and  Kato23  for  both  surfaces,  although  due  to  the  lack  of 
information  it  has  been  assumed  that  there  is  no  angular  dependence  to  the  long 
range  attractive  potential  (i.e.,  the  parameter  y  in  Equation  (6)  is  independent  of 
bond  angle). 

3.  Classical  Trajectory  Methodology 

The  quasiclassical  trajectory  (QCT)  method29  is  used  to  compute  the  thermal 
reaction  rate  constants,  and  the  translational  energy  dependence  of  the  reaction  cross 
sections  and  the  final  vibrational/rotational  state  distributions.  The  validity  of 
classical  mechanics  for  calculating  reaction  rate  constants  and  product  distributions 
has  been  discussed  elsewhere.29  The  classical  approach  is  expected  to  be  reliable  for 
calculating  state-to-state  relaxation  and  reaction  rate  constants  for  which  the 
transition  or  reaction  probability  is  large  (strongly  classically  allowed).  Calculations 
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are  carried  out  separately  for  the  1 2ANand  1 2AOadiabatic  potential  energy  surfaces 
and  then  combined  after  multiplying  the  individual  rate  constants  by  the  temperature 
dependent  electronic  degeneracy  factors  (i.e.,  the  probability  of  a  collision  occurring 
on  a  particular  surface) 


f(T)  =  2<[5  +  3  exp(-22  7.8/T)  +  exp(-326. 6/1 )]  [2  +  2  exp(-l  77.  l/T)\}-‘ . 


A  thorough  discussion  of  the  spin-orbit  effects  and  reactivity  in  the  similar 
0(3P)  +  0H(2n)  system  is  given  by  Graff  and  Wagner.30 

Standard  Monte  Carlo  techniques  are  used  to  compute  the  vibrational 
relaxation  rate  constants  as  a  function  of  temperature  for  specified  initial  vibrational 
states,  v,  where  the  translational  energy,  Ex,  and  the  initial  rotational  states,  j,  are 
chosen  from  Boltzmann  distributions  at  the  given  temperature.  Final  NO(vl\ljl\| 
distributions  are  also  calculated  as  a  function  of  initial  vibrational  state.  The  effect 
of  electronic  angular  momentum  on  the  NO  rotational  state  has  been  ignored,  as  the 
NO  diatomic  properties  appropriate  to  the  NO(X2n1/2)  ground  state  have  been  used. 
The  final  vibrational  and  rotational  quantum  numbers  (v  1^1  j  IS|  are  obtained  from  the 
correspondence  rules  vl^J,  =  (vlMf2  )h  and  jl^J,  =  (jlMf2  )h,  where  vr^J,  and  jl^  are  the 
classical  vibrational  and  rotational  action  variables.  The  final  (vl\ljl\|  distributions  are 
then  obtained  using  the  standard  histogram  method. 

The  classical  equations  of  motion  are  integrated  using  a  variable  step  size 
predictor-corrector  method.31  Although  this  numerical  integration  scheme  has  been 
shown  to  accurately  integrate  long-lived  classical  trajectories  in  various  systems,31-33 
approximately  5%  of  the  trajectories  forming  collision  complexes  did  not  satisfy  time- 
reversal  criterion  even  though  total  energy  and  angular  momentum  were  conserved. 
However,  a  survey  of  these  trajectories  indicated  that  while  the  individual  trajectories 
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were  inaccurate,  the  properties  of  the  ensemble  of  trajectories  were  insensitive  to  the 
accuracy  of  the  numerical  integration  (i.e.,  the  error  parameters).  Thus,  any  error 
in  the  distributions  reported  here  reflect  the  statistical  uncertainty  of  the  calculation. 
A  total  of  38,000  (f9,000  for  each  PES)  trajectories,  of  which  approximately  1/2  (1/4) 
formed  collision  complexes  on  the  2AN(2AQ  PES,  were  used  in  the  current  study. 

4.  Results 

The  room  temperature  0+N0(v)  relaxation  rate  constants  are  compared 
with  the  available  experimental  measurements  in  Figure  3.  While  the  calculated  v=l 
to  vINFO  rate  constant  using  the  two  lowest  (attractive)  potential  energy  surfaces  is 
approximately  a  factor  of  2  below  the  experimental  measurements  of  Fernando  and 
Smith13  and  Lilenfeld,14  the  OCT  deactivation  rate  constants  for  v=2  and  3  are  30- 
40%  greater  than  the  measurements  of  Dodd  et  al.15  Although  it  is  clear  that  the 
present  calculations  are  qualitatively  consistent  with  the  laboratory  data,  it  is 
important  that  the  discrepancy  between  the  v=f 13,14  and  the  v=2,315  measurements 
be  resolved  because  of  the  impact  on  the  modeling  of  NO  emission  in  the 
thermosphere. 

It  is  clear  from  the  present  QCT  calculation  that  the  vibrational  relaxation  on 
both  the  2AN  and  2AO  potential  energy  surfaces  occurs  exclusively  through  the 
formation  of  the  N02  complex.  Therefore,  it  is  interesting  to  compare  the  present 
results  with  the  statistical  adiabatic  channel  model  of  Quack  and  Troe17  which  was 
developed  for  unimolecular  dissociation.  Their  approach  uses  known  molecular 
properties  and  high  pressure  recombination  rate  constants  of  N02  to  calculate  the 
temperature  dependence  of  NO(v)  relaxation  by  O  atoms.  At  300  K,  their  model 
predicts  relaxation  rate  constants  of  1.7x10  11  cm'/s  and  2.0x10  11  cm3/s  for  NO  v=l 
and  v=2,  respectively.  These  values  are  approximately  50%  lower  than  the  QCT 
relaxation  rate  constants  of  2.8x10  11  cm3/s  and  3.1x10  11  cm3/s  for  NO  v=l  and  v=2. 
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However,  Quack  and  Troe  assumed  that  only  the  ground  state  2ANPES  participates 
in  the  relaxation  process.  Based  on  the  recent  work  of  Katagiri  and  Kato,  the  2AO 
PES  should  also  be  considered  in  calculation  of  the  relaxation  rate  constant  which 
would  increase  the  Quack  and  Troe  result  by  up  to  a  factor  of  2.  For  the  QCT 
calculation,  the  2AOPES  contributes  approximately  25%  to  the  relaxation  rate 
constant,  which  indicates  the  present  results  for  NO(v=f,2)  are  in  excellent 
agreement  with  the  statistical  model  of  Quack  and  Troe  for  the  ground  state  surface. 

Figure  4  shows  the  calculated  temperature  dependence  of  the  v=f  to  vN=0 
relaxation  rate  coefficient  compared  to  the  experimental  measurements.  As 
discussed  previously,  the  calculated  room  temperature  results  are  up  to  a  factor  of 
~2  too  low  compared  to  experiment.  However,  the  comparison  of  the  high 
temperature  results  with  the  measurements  of  Glanzer  and  Troe16  are  in  excellent 
agreement.  Although  the  temperature  dependence  of  the  QCT  relaxation  rate 
constant  shows  less  than  a  5%  variation  over  the  temperature  range  of  300  K  to  2700 
K,  the  rate  constant  calculated  without  including  the  electronic  degeneracy  factor 
actually  increases  by  a  factor  of  1.65.  The  effect  of  including  the  temperature 
dependence  of  the  electronic  degeneracy  factor  for  the  2ANand  2AOsurfaces  defined 
in  Equation  (fO),  which  decreases  from  0.095  at  300  K  to  0.060  at  2700  K,  is  to 
produce  a  rate  constant  which  is  independent  of  temperature.  As  temperature 
increases,  more  collisions  occur  on  the  higher  lying  repulsive  potential  energy 
surfaces  which  are  ineffective  in  vibrational  relaxation.  Thus,  it  is  important  to 
distinguish  between  the  temperature-dependent  effects  due  to  the  population  of 
electronic  states  and  the  vibrational  relaxation  on  the  particular  PES.  Clearly,  for 
systems  as  complicated  as  O+NO,  establishing  a  mechanism  (i.e.,  statistical  behavior) 
for  vibrational  relaxation  can  not  be  decided  by  the  overall  temperature  dependence 
of  the  rate  constant  alone.  It  is  interesting  to  note  that  the  relaxation  rate  constant 
predicted  by  the  statistical  adiabatic  channel  model17  decreases  from  f. 7x10  11  cm3/s 
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at  300  K  to  1.2x10  11  cm3/s  at  2100  K  in  disagreement  with  the  OCT  temperature 
independent  relaxation  rate  indicating  that  the  present  calculations  do  indeed  have 
a  significant  nonstatistical  component. 

There  is  currently  only  limited  experimental24  and  theoretical  information17 
concerning  the  importance  of  multiquantum  versus  single  quantum  relaxation  for 
0+NO(v).  The  ab  initio  electronic  structure  calculations,23  which  indicate  that  the 
NO  vibrational  relaxation  occur  on  surfaces  dominated  by  a  strongly  bound  complex, 
would  argue  for  efficient  multiquantum  relaxation. 17,32  In  fact,  it  has  been  shown  by 
Glanzer  and  Troe  for  NO(v=2)  at  high  temperatures  (2700  K)  and  Fernando  and 
Smith  for  NO(v=l)  at  room  temperature  that  the  vibrational  deactivation  is 
consistent  with  the  formation  of  a  collision  complex.  In  addition,  the  recent 
measurements  of  Dodd  et  al.24  for  v=3  to  vlNF2  relaxation  rate  coefficient  confirms 
the  expectation  of  multiquantum  relaxation.  However,  no  information  on  the 
relaxation  of  highly  vibrationally  excited  NO  currently  exists.  Detailed  QCT 
calculations  of  the  vibrational  relaxation  of  NO(vo  9)  by  oxygen  atoms  at  300  K  are 
shown  in  Figure  5.  The  results,  which  indicate  that  multiquantum  relaxation  is  very 
efficient,  are  independent  of  final  state  and  can  be  qualitatively  represented  by 
kYV- =  kv/v  ■  Furthermore,  the  present  results  are  in  qualitative  agreement  with  the 
results  of  the  statistical  adiabatic  channel  model17  for  NO(v=l,2)  and  the 
expectations  of  statistical  theories17,33,34  which  predict  that  vN=0  should  be  the  most 
populated  final  vibrational  state  independent  of  initial  v.  Finally,  the  experimental 
branching  fraction  of  0.35"  0.10  for  v=3  to  vN=2  measured  by  Dodd  et  al.  is  in 
excellent  agreement  with  the  QCT  calculated  value  of  0.29"  0.02  providing  further 
evidence  of  efficient  multiquantum  relaxation  in  O+NO  collisions. 

It  is  difficult  to  establish  whether  a  particular  dynamical  system  can  be 
described  by  statistical  theories.  Wagner  and  Parks34  have  discussed  in  detail  a 
general  framework  encompassing  classical  statistical  theories,  with  emphasis  on  the 
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concept  of  a  strong  coupling  region.  This  approach  was  employed  in  previous  work33 
investigating  the  dynamical  origin  of  statistical  behavior  in  classical  systems  which 
exploited  the  sensitivity  of  classical  dynamics  to  initial  conditions.  It  was  argued  that 
exponential  divergence  of  trajectories  initially  adjacent  in  phase  space  was 
responsible  for  memory  loss  in  bimolecular  collisions.  Since  such  a  detailed  study  is 
beyond  the  scope  of  the  present  work,  one  can  use  simple  criterion  for  statistical 
behavior  in  the  O+NO  relaxation,  for  example,  the  ratio  of  the  reactive  (i.e.,  O  atom 
exchange)  to  nonreactive  relaxation  rate  constant.  In  a  symmetric  system,  such  as 
O+NO  where  the  reactants  and  products  are  identical,  any  statistical  theory  would 
predict  the  ratio  of  the  nonreactive  to  reactive  relaxation  rate  constant  to  be  unity.33 
This  expectation  implicitly  assumes  that  the  potential  energy  surface  is  symmetric 
with  respect  to  the  interchange  of  reactants  and  products.  The  simple  test  of  this 
hypothesis  for  the  0+N0(v)  relaxation  is  illustrated  in  Figure  6,  where  the  reactive 
and  nonreactive  components  of  the  vibrational  relaxation  rate  constants  presented 
in  Figure  3  are  shown.  The  OCT  results  indicate  that  the  nonreactive  relaxation 
process  is  approximately  a  factor  of  2  more  efficient  than  the  reactive  O  atom 
exchange  relaxation  for  all  initial  vibrational  states.  This  result  implies  that  although 
the  NO(v)  relaxation  proceeds  via  a  collision  complex  and  vibrational  energy  may  or 
may  not  be  equilibrated,  the  large  fraction  of  the  collisions  which  favor  reforming 
reactants  indicates  significant  nonstatistical  behavior.  Similar  results  have  been 
observed  in  a  molecular  beam  study35  of  the  reactions  of  alkali  atoms  with  alkali 
halides  and  in  a  QCT  study34  of  alkali  halide  exchange  reactions.  Thus,  the  extent  of 
statistical  behavior  in  0+N0(v)  relaxation  remains  to  be  firmly  established. 

5.  Summary 

Detailed  quasiclassical  trajectory  calculations  of  the  0+N0(v)  vibrational 
relaxation  rate  constants  have  been  carried  out  using  the  two  lowest  electronic 
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potential  energy  surfaces  neglecting  nonadiabatic  interactions.  The  calculated  total 
relaxation  rate  constants  are  in  good  agreement  with  the  available  experimental  data 
and  the  statistical  adiabatic  channel  model  of  Quack  and  Troe.  Furthermore,  the 
QCT  calculations  of  the  state-to-state  relaxation  rate  constant  show  the  importance 
of  multiquantum  vibrational  energy  transfer,  again  consistent  with  statistical 
behavior,  and  in  excellent  agreement  with  the  recent  experimental  measurement  of 
the  v=3  to  vl\P2  branching  ratio  by  Dodd  et  al.  Although  the  relaxation  is  shown  to 
proceed  exclusively  through  the  formation  of  a  long-lived  bound  intermediate,  and 
the  state-to-state  rate  constants  are  in  good  agreement  with  expectations  based  upon 
statistical  theories,  the  branching  ratio  (i.e.,  nonreactive  vs.  O  atom  exchange) 
indicates  that  the  collision  complex  retains  some  memory  of  the  initial  state. 

The  present  results  and  the  recent  measurements  of  Dodd  et  al.  for 
NO(v=2,3)  imply  that  the  0+N0(v=f)  relaxation  rate  constant  may  be  a  factor  of 
2-3  smaller  than  previously  thought.  Reducing  the  NO(v=f)  relaxation  rate  constant 
used  in  NO  thermospheric  cooling  rate  calculations  would  require  higher 
atmospheric  temperatures  and/or  an  increase  in  the  NO/O  number  densities  in  the 
models  to  maintain  agreement  with  measurements.  The  suggestion  of  higher 
thermospheric  NO  densities  for  cooling  rate  calculations  is  consistent  with  numerous 
measurements  of  NO  densities  in  the  lower  thermosphere  which  show  that  measured 
values  are  generally  larger  than  those  predicted  by  photochemical  models. 
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TABLE  1.  Parameters  for  the  LEPS  Potential1 


De 

re 

ai 

a2 

a3 

z 

NO(x2n1/2) 

6.6144 

1.1508 

5.035 

5.151 

2.998 

0.3 

02(X2X+g) 

5.2132 

1.2075 

5.476 

7.950 

6.341 

0.2 

aThe  diatomic  parameters  are  taken  from  ref.  25.  Energies  are  given  in  eV  and  distances  in 
Dngstroms. 


TABLE  2.  Parameters  for  the  N02  f  2ANand  f  2AOExponential  Functions3 


1 2AN  1 2  AO 


c0 

19.00085 

48.52510 

Cl 

-62.34160 

-194.92802 

c2 

52.40209 

245.46673 

C3 

-9.43063 

-93.79654 

y 

2.30 

2.30 

r° 

1.195 

1.270 

“Energies  are  given  in  eV  and  distances  in  Dngstroms. 
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x  (A) 

Figure  1.  Polar  contour  plot  of  the  ONO  2ANpotential  energy  surface  with  the  NO 
internuclear  distance  set  equal  to  its  equilibrium  value.  The  contour  values  in  eV  are 
4(f),  3(2),  2(3),  1.5(4),  1(5),  0.5(6),  0.25(7),  -0.1(8),  -0.25(9),  -0.5(10),  -1(11),  and  - 
2(12)  relative  to  the  0(3P)+N0(2n1/2)  asymptote. 


x  (A) 

Figure  2,  Polar  contour  plot  of  the  ONO  2AOpotential  energy  surface  with  the  NO 
internuclear  distance  set  equal  to  its  equilibrium  value.  The  contour  values  in  eV  are 
4(1),  3(2),  2(3),  1.5(4),  1(5),  0.5(6),  0.25(7),  -0.1(8),  -0.25(9),  and  -0.5(10)  relative  to 
the  0(3P)+N0(2n1/2)  asymptote. 


18 


Figure  3.  Room  temperature  (300  K)  rate  coefficient  for  the  deactivation  of 
vibrationally  excited  NO(v)  upon  impacts  with  atomic  oxygen  as  function  of  the 
vibrational  level.  The  present  QCT(nMn)  calculation  is  compared  to  the 
experimental  measurements  of  Fernando  and  Smith13  (#),  Lilenfeld14  (cf),  and  Dodd 
et  al.15  (>).  The  rate  coefficient  kv  is  the  sum  of  the  rate  coefficients  for  deactivation 
to  all  the  final  levels,  kv  =  ^  kw •  ■ 
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Figure  4.  Temperature  dependence  of  the  vibrational  deactivation  of  NO(v=l)  by 
atomic  oxygen.  The  present  QCT(nMn)  calculation  is  compared  to  the  experimental 
measurements  of  Fernando  and  Smith13  (#),  Lilenfeld14  (et),  and  Glanzer  and  Troe16 
(>)• 
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Figure  5.  Rate  coefficients  for  the  deactivation  of  NO(v)  to  NO(vl\|  by  atomic  O  as 
a  function  of  the  final  vibrational  level  vN  The  final  vibrational  levels  are  connected 
by  solid  lines  for  clarity  of  presentation  and  the  initial  level  v  is  indicated  on  each 
curve. 


Figure  6.  The  contribution  of  the  nonreactive  (n"  n)  and  O  atom  exchange  (riMn) 
reactions  to  the  rate  coefficient  for  vibrational  relaxation  kv  at  room  temperature  as 
a  function  of  the  initial  vibrational  level  v. 
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